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Now you see it, now you don't - the circumstellar disk in 
the GRO J1008-57 system 
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ABSTRACT 

Multiwavelength observations are reported here of the Be/X-ray binary pulsar 
system GRO J1008-57. Over ten years worth of data are gathered together to show 
that the periodic X-ray outbursts are dependant on both the binary motion and the 
size of the circumstellar disk. In the first instance an accurate orbital solution is 
determined from pulse periods, and in the second case the strength and shape of the 
Ha emission line is shown to be a valuable indicator of disk size and its behaviour. 
Furthermore, the shape of the emission line permits a direct determination of the 
disk size which is in good agreement with theoretical estimates. A detailed study 
of the pulse period variations during outbursts determined the binary period to be 
247.8±0.4d, in good agreement with the period determined from the recurrence of the 
outbursts. 
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• H , 1 INTRODUCTION AND BACKGROUND 

5^ 1 The Be/X-ray systems represent the largest sub-class of 
' massive X-ray binaries. A survey of the literature reveals 
that of the 115 identified massive X-ray binary pulsar sys- 
tems (identified here means exhibiting a coherent X-ray 
pulse period), most of the systems fall within this Be coun- 
terpart class of binary. The orbit of the Be star and the com- 
pact object, presumably a neutron star, is generally wide and 
eccentric. X-ray outbursts are normally associated with the 
passage of the neutron star close to the circumstellar disk 
(Okazaki & Negueruela 2001). A review of such systems may 
be found in Coe et al. (2000). 

The source that is the subject of this paper, the X-ray 
transient GRO J 1008-57, was discovered on July 14 1993 
by the BATSE experiment on the Compton Gamma Ray 
Observatory (Stollberg et al. 1993 and Wilson et al. 1994). 
The source was observed in the 20-200 keV band, and was 
found to be pulsating at a period of 93.5s. Its spectrum 
was observed to be consistent with optically thin thermal 
bremsstrahlung with kT=25 keV. From these X-ray data it 
was concluded that the system was a massive binary system - 
a neutron star with either a Be or supergiant primary. A full 
report of this discovery outburst was presented in Shrader 
et al. (1999). 



Subsequently the optical counterpart was identified by 
Coe et al, (1994) who showed it to be a V=15.3 OB star 
with a strong IR excess and strong Balmer line emission. In 
this paper we present more than 10 years of optical monitor- 
ing of this counterpart, including a blue-end spectrum taken 
from the Southern African Large Telescope (SALT), which 
enables refinement of the spectral class to much higher ac- 
curacy than previous published work. We also demonstrate 
how the optical characteristics correlate with the long-term 
X-ray (RXTE/ASM data) behaviour. This is supplemented 
by 7-ray observations from INTEGRAL, in particular, re- 
porting the details of the June 2004 outburst from GRO 
J 1008-57. 



2 OPTICAL DATA 
2.1 Red spectra 

Ha data have been collected over the last 13 years from 
a series of telescopes. The dates and properties of the Ha 
line are presented in Table [1] In this table the following 
telescopes and configurations have been used: 

• AAT - 3.9m telescope, AAT (Australia), Royal Green- 
which Observatory Spectrograph (RGOS), 25cm camera, 
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Table 1. Tabic of Ha measurements. See text for details of ob- 
servatory/instrument used. 



Observation 


Observatory 


Ha EW 


Peak 


Date 




(A) 


Separation 

(km s~ ) 


21 Dec 1993 


SAAO 


-17±2 




24 Dec 1993 


SAAO 


-13.6±1.3 


250 


26 Feb 1994 


AAT 


-19.4±1.2 


215 


07 Mar 1994 


SAAO 


-19.6±0.5 




26 Feb 1995 


SAAO 


-19±3 




19 Nov 1995 


SAAO 


-11±2 




20 Jun 1997 


SAAO 


-2.6±1.5 




20 Jun 1997 


SAAO 


-4.9±1.3 


350 


06 Feb 1998 


SAAO 


-6.7±1.2 




11 Jan 1999 


SAAO 


-19.0±1.6 




14 Dec 2002 


SAAO 


-20±2 




15 Mar 2005 


SAAO 


-22.6±1.7 




19 Mar 2005 


SAAO 


-21.5±1.5 


243 


07 May 2006 


SAAO 


-20.9±1.2 





1200V grating and TEK (1024x1024) CCD. The dispersion 
was 0.8A/pixel and the signal to noise ratio from the 600s 
exposure was ~40. 

• SAAO - 1.9m telescope, Sutherland Observatory (South 
Africa), spectrograph, SITe detector, 1200 1/mm grating. 
The dispersion was l.oA/pixel and the signal to noise ra- 
tio ~10. 

Selected Ha line profiles are presented in Figure [1] 
where more than one spectra exists for a particular year 
then an average of the two is shown. A clear evolution of 
the Balmer line over the 10 year period is seen. Apart from 
the structural changes which will be discussed below, the 
strength of the Ha emission line varies considerably over 
the ~10 year period - see Figure [2] . The equivalent width 
determinations can be seen to go through a significant de- 
cline in the middle of the period, recovering strongly in the 
latest data. 

The peak separations were determined by fitting 
Lorentzian profiles to the Ha emission lines. In the cases 
where the x 2 was significantly improved by including two 
lines, as opposed to one line, the line parameters were used 
to measure the separation (in velocity space) between the 
red and violet peaks. These kind of changes in the line pro- 
files are common in Be stars (see, for example, Telting & 
Kaper 1994). 

2.2 Blue spectrum 

A blue end spectrum was obtained from SALT on 18 June 
2006 - see Figure [3] The spectrum was obtained using the 
Robert Stobie Spectrograph in long slit mode, with a grating 
resolution of 1800 1/mm and an exposure time of 600s. This 
gave a dispersion of 0.34A/pixel and a signal to noise varying 
from ~10 (blue end) to ~50 (red end). 

Mosaicking of the three chips was performed by using 
the SALT IRAF packages (the two chip gaps are clearly 
visible in Figure [3} as part of the pipeline processing of the 
data. Further wavelength calibrations were performed using 
an Ar arc lamp solution. No flux calibration has been made. 

The prominent feature is the strong double peaked H/9 
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Figure 1. The evolution of the Ha profile over 13 years. Individ- 
ual spectra have been arbitrarily re-scaled in flux to permit the 
shape of the Ha profile to be clearly visible. 



line, which has an equivalent width of —1.74 ± 0.09A. The 
only previous blue end spectrum (Coe et al.1994) gives a H/3 
equivalent width of — 1.4±0.lA. Two other Balmer lines (H7 
and He) are visible; but H<5 seems to have been effectively 
lost due to infilling. This is unusual given that the Balmer 
lines on either side do not seem to be so affected. The spec- 
trum has been carefully checked for any evidence of scat- 
tered light or other instrument atrtefacts, but no evidence 
for any such features were found. The lack of any identifi- 
able Hell lines, in particular 4686 A indicates that the star 
is Bl or later (Walborn and Fitzpatrick 1990). Furthermore 
the presence of the Oil CIII blend at 4650A restricts the 
classification to the range B1-B2 (Walborn and Fitzpatrick 
1990). It is worth noting, however, that the lack of other 
metal lines is a little surprising but may be due, in part, to 
the decreasing signal-to-noise ratio towards the blue end of 
the spectrum. However, such a lack of lines is not uncommon 
in low metalicity environments such as the Small Magellanic 
Cloud (McBride et al, 2007). 



3 ORBITAL X-RAY PROFILE 

X-ray monitoring data were collected by the All Sky Moni- 
tor (ASM) instrument on the Rossi X-ray Timing Explorer 
(RXTE) observatory. The data consist of daily flux aver- 
ages covering the X-ray energy range 1.3 - 12.1 keV. There 
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Figure 2. The evolution of the flux of the He* emission line over 
13 years (1993-2006). 



1000 




4000 



4200 4400 
Wavelength (A) 



4600 



4800 



Figure 3. SALT blue spectrum of the optical counterpart to 
GRO J1008— 57. Major features in the spectrum are indicated. 
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Figure 4. BATSE (top panel) and RXTE/ASM (lower two pan- 
els) data folded at the period of 248.9d (Levine & Corbet, 2006) 
for three different time intervals. The phases shown are those cor- 
responding to the ephemeris given in Section 6.1. In each case a 
±lcr error bar is shown. 



at the same period and same phase. These data cover the 
harder spectral range of 20 - 100 keV, but reveal a very 
similar modulation occurring in the period 1991-1996 - see 
upper panel in Figure [4] The phase of the outburst and the 
degree of modulation looks very similar to that seen with 
RXTE/ASM in recent years. 



is no evidence in the daily data averages for any significant 
detection of the source on short timescales (~ a few days). 

The ~10 years worth of data were folded at the pro- 
posed binary period of 248. 9d determined by Levine & Cor- 
bet, 2006 from the same data. A strong, clear modulation 
at that period is very obvious. 

Since it is clear from Figure [2] that the strength of the 
Ha line has been significantly changing over the 10 year 
interval of the ASM data, the complete data set was divided 
into two equal periods (1996-2001 & 2001-2006). Each data 
set was separately folded at this same period and with the 
same phase. The results are shown in the lower two panels of 
Figure [4] from which it is clear that the degree of modulation 
had changed very significantly over the 10 years. 

To explore the modulation at even earlier times (i.e. 
before 1996) data from the BATSE occultation instrument 
on Compton Gamma Ray Observatory (CGRO) were folded 



4 CGRO /BATSE ORBITAL SOLUTION 

The data types that are available from the 8 BATSE detec- 
tors on CGRO are the 16-channel continuous or "CONT" 
data, sampled at 2.048s intervals, and the 4-channel dis- 
criminator or "DISCLA" data sampled every 1.024s. For 
the analyses in the present work the DISCLA data were used 
from which pulsed fluxes and frequencies were extracted (see 
Bildsten et al.1997 for a full explanation on the use of DIS- 
CLA data for pulsar observations). 

To assess the confidence in the detection of a known pul- 
sar, the data were searched through in a narrow frequency 
range centred on the pulsar's nominal frequency. 

First the DISCLA rates were grouped into 300 s inter- 
vals. A pulse profile was then estimated using a 6th order 
Fourier expansion for each of these intervals. These profiles 
were then grouped into 4-day intervals. Within each 4-day 
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interval, the Fourier coefficients for each 300-s pulse profile 
were re-computed for a grid of offsets from the nominal fre- 
quency. The best frequency offset was determined using the 
Z„ statistic (Buccheri 1983), expressed as: 



E 



(i) 



are the mean Fourier coefficients and <r Mj! their 



where [ik 

Poisson errors; TV can take any value from 1 to 6 (in general 
TV = 3 or 6 was used). In addition, within each 4-day interval, 
a mean pulse profile was determined for each frequency off- 
set by chi-square minimization. Unfortunately, BATSE noise 
is generally non-Poissonian (Bildsten et al. 1997), so <r M( . 
is multiplied by the reduced chi-square factor to produce 
the new statistic Y n . Monte-Carlo methods are required to 
convert Y n into percentile significances in order to establish 
confidence levels for the detected signal. An advantage to 
using this method is that any number of harmonics can be 
excluded from the statistic if there are other sources within 
the field of view with pulsations at those particular harmonic 
frequencies. 

If an outburst lasts for a significant fraction of a sys- 
tem's binary orbital period, the size and shape of the orbit 
can be determined by analysing Doppler shifts and pulse ar- 
rival delays in the pulsar signal. As an outburst progresses, 
BATSE obtains a continuous history of the daily source in- 
tensity and frequency; this detected frequency will be the 
intrinsic spin frequency of the pulsar plus or minus a shift 
due to the neutron star's orbit around its companion. The 
process of calculating the orbit is further complicated by the 
fact that the neutron star is usually spinning up or down 
during an outburst, and so it is necessary to dissentangle 
the orbital from the spin torque effects. 

The method used with the BATSE data is the piece- 
wise approach explained in detail in Wilson et al, 2003. In 
essence, it assumes that the torque during an outburst will 
be approximately constant during small time intervals, such 
that the intrinsic spin up will be linear during this time. 
As such, the data available during an outburst will be di- 
vided into short segments of 3-5 data points each. In order 
to compute the orbital parameters of the system, the data 
in all the segments are fitted with global values of P or t, t p , 
e, ui and a x sini, but a different Vintr in each segment. 

GRO J1008-57 was first detected during a bright 
outburst in July 1993 by BATSE with pulsations at 
93.587±0.005s (Stollberg et al, 1993). This outburst can be 
seen around JD 2449200 in the second from top panel in 
Figure [5] After this initial outburst, a number of less sig- 
nificant detections were made, which could correspond to 5 
further outbursts. These are spaced ~260d apart, which led 
to this being proposed as the orbital period of the system 
(Finger et al.1994). Shrader et al. (1999) find weak evidence 
for a ~135d period in the ASM data from RXTE. 

Orbit fitting calculations were carried out using the 
>99.9% significance data points (all the points illustrated in 
Figure(5]l, producing an orbital period of 247.8±0.4d. This is 
in good agreement with the Levine & Corbet (2006) value of 
248.9±0.5d. The results can be seen graphically in Figure [5] 
and the orbital parameters are presented in Table [2] 

Using the frequencies corrected for orbital motion from 
Figure [5] we derive for the initial outburst a P = —5.6 x 
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Figure 5. BATSE data of GRO J1008-57. Top panel shows 4-day 
averaged frequencies and the second panel down shows the pulsed 
flux (n.b. only the Y n >99.9% significance data points are shown). 
Third panel: the crosses are the emitted frequency returned by 
fitting the model described in the text; the solid line illustrates the 
magnitude of the orbital effect and the circles are the frequencies 
as measured by BATSE. Lowest panel : pulse frequency residuals 
for the model. 



Table 2. Orbital parameters for GRO J1008-57. 



Porbital (d) 

e 

a x sini (light-s) 

Tperiastron (MJD) 



247.8±0.4 
0.68±0.02 
530±60 
49189. 8±0.5 



lO^ss -1 implying L x = 1.0 x 10 36 erg s" 1 (B = 1.4xl0 13 G). 
Compare to L x = 4.1 x 10 36 ergs _1 (B = 6.0 x 10 13 G), ob- 
tained by Shrader et al. 1999. 



5 INTEGRAL 7- RAY DATA 

GRO J1008-57 was detected by INTEGRAL on two occa- 
sions - June 2004 and October 2005 (Grebenev et al., 2005 
have reported on the latter outburst). On many other occa- 
sions the source has fallen within the INTEGRAL field of 
view. The data from the IBIS imaging instrument on-board 
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Folded period: 21504960.0000000 s 



1.317x10* 1.318x10* 1.319x10* 

Time (JD-2440000.5) 



Figure 7. INTEGRAL IBIS data showing the outburst of June 
2004 in the energy range 18—60 keV. 



Start Time 52706 1 8:02:1 4:784 Stop Time 53826 1 1 :36:50:784 



Figure 6. All available INTEGRAL data in 18-60 keV range 
folded at the outburst period of 248. 9d. The phase is the same as 
in Figure [4] 



INTEGRAL were reduced and analysed using the standard 
Offline Science Analysis, OSA, software version 5.1. 

A light curve was obtained by analysing all available IN- 
TEGRAL data when the source was within the IBIS/ISGRI 
field of view (1085 science windows) and extracting the 18-60 
keV flux at the source position in each science window. The 
sampling time for this light curve is therefore at the typical 
science window duration of ~2000s. Folding the light curve 
on the ephemeris given in Section 6.1 yields the folded light 
curve shown in Figure [6] Though the modulation shown is 
dominated by the one major outburst of June 2004 the pro- 
file is in good phase agreement with that shown in Figure 2] 
for the BATSE and RXTE/ASM outbursts. 

Since there are not many data on the second outburst, 
the results from the first, previously unreported outburst are 
presented in this section. 

A lightcurve of the 2004 outburst is shown in Figure [7] 
for the photon energy range 18-60 keV. From this figure it 
is possible to ascertain that the main burst lasted less than 
23.5 days, and rose in an approximately linear manner with 
a doubling time of ~3.5d. In total INTEGRAL observed this 
outburst for ~7.9 d. 

A 10s binned light curve in the 20-40 keV photon en- 
ergy range was obtained from the IBIS /ISGRI data using the 
OSA iijight tool. The finely binned lightcurve was analysed 
using the Lomb-Scargle periodogram method by means of 
the fast implementation of Press & Rybicki (1989). The re- 
sulting power spectrum shows a clear strong peak at 0.01068 
Hz with a power of ~435. Using the error formulae derived 
by Home & Baliunas (1986) we find that this corresponds to 
a periodicity of 93. 6672T0. 0005s. This period was indepen- 
dently confirmed by the 4-10 keV data from the X-ray cam- 
era on-board INTEGRAL, JEM-X. Furthermore, the period 
is in very close agreement to that of 93.62T0.01 reported by 
Shrader et al. (1999) from the 1993 ASCA data, indicating 
very little spin up (or down) over more than a decade. 

The phase folded pulse profile was determined in 4 sep- 




Figure 8. INTEGRAL IBIS and JEM-X pulse profiles. 



arate energy channels; 4-10 and 10-20 keV from JEM-X 
and 20-40 and 40-60 keV from IBIS/ISGRI. The resulting 
single peaked profiles are shown in Figure [5] From these 
profiles it was possible to determine the pulse fractions, 
(Imax — Imin)/I ma x, to be 72±8% in the 4-10 keV band, 
70±10% in the 10-20 keV band, 80±10% in the 40-60 keV, 
and 100±30% in the 40-60 keV range. The IBIS/ISGRI 
pulse fractions must be treated with care as they are derived 
from the folded fine timing lightcurve and consequently are 
not absolutely flux calibrated. 

Spectra for the 2004 outburst of GRO J1008-57 were 
produced using OSA for both the IBIS/ISGRI and JEM-X 
instruments. The source was not always within the field of 
JEM-X and consequently there were only 12 science windows 
of data available spread across the outburst where the source 
was detected by that instrument. The IBIS/ISGRI data for 
the same time span comprised of ~80 science windows. The 
spectra were fit simultaneously in XSPEC 12.3.0; the JEM- 
X data were fit over the 3-20 keV range and the IBIS/ISGRI 
data over the 20-100 keV range. A cross-calibration con- 
stant, C isgri, was added to all spectral models to account 
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Table 3. The results of the simultaneous spectral fits of the JEM- 
X and ISGRI data. 



Parameter 



Cut-off Power Law Thermal Bremsstrahlung 



C ISG ri 1-6 ± 0.1 

cut-off energy 8.0 ± 1.0 keV 

Fold energy 21 ± 2 keV 

Photon index 1.4 ± 0.1 
kT 

X 2 /dof 150.46/146 



1.5 ± 0.1 



21 ± 1 keV 
162.31/148 




be defined as phase 0.0. From this the following outburst 
ephemeris is determined: 

T outburst = 50186±4 + 248.9n 
where the dates are in Modified Julian Day format and n is 
the outburst number. 

If the time of outburst is compared to the predicted 
time of periastron determined from the pulse period varia- 
tions (see Table 2), it agrees within 1.5d for the MJD 49189 
outburst, indicating no significant phase shift between these 
two events. 

If we use this RXTE/ASM ephemeris then we can de- 
termine that the phase of the two INTEGRAL outbursts 
discussed here are 0.99 and 0.00. Furthermore, the phase of 
the peak of the 1993 BATSE outburst reported in Shrader 
et al. (1999) is 0.01. Thus all three reported events fit very 
comfortably with the above ephemeris and hence are almost 
certainly Type I outbursts. This is supported by comparing 
the phase of the peaks of outbursts as seen in the folded 
lightcurves from BATSE, RXTE/ASM and INTEGRAL all 
of which seem to match very well. There is no evidence for 
the phase changes in outburst seen in EXO 2030+375 and 
explained by Wilson et al. (2002) as related to global one 
arm oscillations in the circumstellar disk. Though there are 
changes in the disk, evidenced by the overall shape changes 
of the Ha profile seen in Figure [1] they are obviously not 
great enough to result in such outburst phase shifting. 



Figure 9. INTEGRAL combined IBIS and JEM-X spectrum ob- 
tained by deconvolving the data using the Cut-off Power Law 
model. 



for any uncertainties in the matching of the two instruments. 
A systematic error of 2% was used. 

Two spectral models were fit to the data, a power-law 
with a high energy cut-off and a thermal bremsstrahlung 
model. Both models fit equally well and have a reduced chi- 
square of ~1; the details of the fit parameters are shown 
in Table [3] An attempt was made to introduce an absorp- 
tion component to both, however the fit could not resonably 
constrain the parameter. Consequently, an absorption com- 
ponent fixed to the anticipated galactic absorbing column 
density of ~1.5 x 10 22 cm -2 was added. The unfolded spec- 
tra with the high-energy cut-off power law model fitted is 
shown in Figure [9] 

The fluxes given by the spectral models were 5.7 x 
10~ 10 erg cm -2 s _1 in the 3-20 keV band and 4.1 x 10~ 10 
erg cm -2 s _1 in the 20-100 keV band. Taking these numbers 
combined with the source distance of 5kpc (Coe et al, 1994) 
produces a flux estimate of (1-2) x 10 36 erg s _1 . This value 
is in good agreement with the level of X-ray luminosity es- 
timated from the small accretion-torque driven pulse period 
changes reported in Section 4. 



6 DISCUSSION 

6.1 X-ray outburst ephemeris 

Using the clear outburst profile presented in Figure [4] for 
the period 2001-2006 it is possible to determine an accurate 
ephemeris. We assume the binary period is 248. 9d (Levine 
& Corbet, 2006) and chose the highest bin in the figure to 



6.2 Determining the circumstellar disk size and 
its relationship to X-ray activity. 

Figure [5] shows the Ha equivalent width measurements over 
the last 13 years. Although the Ha equivalent width is not a 
direct measurement of the size of the circumstellar disk the 
data in this figure can be used as an indicator of growth and 
decline in the disk. To this effect we interpret the significant 
decline in the Ha flux, followed by a strong recovery as a 
period of disc shrinking and regrowth. This idea is consistent 
with the lack of observed X-ray outbursts of GRO J1008— 57 
during 1996-2001 - see Figured 

Another feature to take note of in Figure [5] is the maxi- 
mum measured Ha equivalent width. Reig, Fabregat & Coe 
(1997) show that this maximum Ha equivalent width is cor- 
related with the orbital period in a Be/X-ray binary and the 
observed correlation can be explained by a viscous circum- 
stellar disk truncated by the resonant torque of the orbiting 
neutron star (Negueruela et al., 2001). The system under dis- 
cussion here spends a large fraction of the observed timescale 
at an Ha equivalent width of close to — 20A, while under- 
going Type I X-ray outbursts. This may suggest that the 
—20 A EW measurement describes the circumstellar disk at 
a size that puts it close to the L\ point at periastron. 

To estimate the size of the Ha emitting region of the 
circumstellar disk we make used of the peak separation mea- 
surements in Table 1 and the relation from Huang (1972). 

where r is the disk radius, i is the inclination and AV 
is the peak separation. The average of the three measure- 
ments taken while the system was in normal outburst was 
determined to be 236±15 km/s and used in the following 
determination of the disk size. 

From the blue spectrum we determine the spectral type 
of the optical counterpart of GRO J1008-57 to be Bl - B2 
Ve. From Allen (1973) we estimate the mass and radius of 
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such a star to be 15 Mq and 7R© respectively. Then, using 
the measured values of a x sin i and assuming that the Be 
star disk is in the plane of the orbit we find an inclination 
angle of i — 36° for the circumstellar disk of GRO J1008— 57. 
These stellar parameters give an estimated Ha disk size of 
72 Rq and a ratio of disk size to Li point (measured from 
the centre of the Be star) of Rj/Rl ~ 0.64 

In comparison, if the data for 1997 are used with a AV 
of 350 km/s, then this results in a much smaller disk of the 
order of 33 Rq in size. 

Even the 72 R© is much smaller than the truncation 
radius predicted by the viscous disk theory for this Be/X- 
ray binary by Okazaki & Negueruela (2001). Using slightly 
different stellar and orbital parameters those authors predict 
that the circumstellar disk of GRO J1008-57 is most likely 
truncated at the 7:1 or 8:1 resonance radius. At such a radius 
the disk verges on the Li point at periastron, thus explaining 
the frequently observed Type I outbursts. 

Exploring the upper and lower limits of the disk size 
estimates presented in this work requires the identification of 
the major uncertainties affecting the results. The uncertainty 
in the spectral types introduces a range of possible masses 
for the Be star of 10 - 18 Mq. Using extreme upper limits 
for both the stellar mass and a^sini, with the lower limit 
on AV, we find that Rd/Rz, ~ 0.86. At the other end of the 
range we find R d /Ri ~ 0.51 

However, it is worth bearing in mind that, contrary to 
the assumption in most disk models, the Be star disk is prob- 
ably not isothermal. Thus the physical size of the disk is not 
constrained to the Ha emitted size of the disk. For example, 
in the case of the isolated Be star Per the Ha disk size 
has been measured by direct interferometry to be around 
ten stellar radii (Quirrenbach et al., 1997), whereas similar 
radio observations (Dougherty & Taylor 1992) suggest that 
material may be present as far away as hundreds and even 
thousands of stellar radii. So, in the case of GRO J1008-57 
the totality of the disk may extend beyond the radius mea- 
sured here from the Ha data. This is indirectly confirmed by 
the very fact that we see regular X-ray outbursts, suggesting 
the disk material must, indeed, reach the Li point. 



7 CONCLUSIONS 

The study of pulse period changes in GRO J1008-57 from 
the BATSE data have allowed a precise orbital determina- 
tion. The orbital period found agrees, within errors, with 
that determined by Levine & Corbet (2006) from X-ray out- 
burst cycles. These outbursts seem to be very reliable in na- 
ture, showing no variation with X-ray energy, nor are they 
correlated with the changes in the Ha profile, as was seen 
in the case of EXO 2030+375 (Wilson et al, 2002). How- 
ever, the substantial changes in the Ha EW correlate nicely 
with the existence and changing levels of the X-ray outburst 
sizes. In addition, using the peak separation measured from 
the Ha profiles an indication was provided of the disk size 
which is in general agreement with the detailed models of 
Okazaki & Negueruela (2001). The results are a direct con- 
firmation of the value of multiwavelength observations of 
these complex systems. 
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